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SUMMARY 


The  theoretical  and  experimental  aspects  of  the  project  have  been 
closely  coordinated  in  keeping  with  the  underlying  philosophy  of  this  project. 
Substantial  progress  has  been  made  theoretically  on  the  development  of  the 
Kramers-Kronig  (Hilbert)  transforms  and  their  application  to  optical  data. 
Kramers-Kronig  analysis  of  optical  reflectivity  data  has  been  used  extensively 
in  the  literature  although  the  reliability  of  this  analysis  until  now  has  been 
questionable.  We  are  now  in  a  position  to  analyze  optical  reflectivity  data 
on  single  crystals  and  thin  films^and  determine  the  fundamental  optical 
constants,  n  and  k,  with  a  reliability  and  accuracy  well  beyond  that  available 
in  the  past. 

In  addition,  computer  programs  are  under  development  which  will  delineate 
the  optimum  experimental  conditions  and  will  perform  analysis  of  the  optical 
constants  of  thin  films  from  transmission  data. 

Basic  to  this  program  is  a  comparison  of  the  properties  of  materials  in 
single  crystal  and  amorphous  form.  The  growth  of  large  single  crystals  of 
several  compounds  and  their  solid  solutions  has  been  successfully  accomplished 
and  a  considerable  amount  of  optical  data  has  been  amassed  on  these  crystals. 

NA  very  significant  finding  ,has  £een  that  the  standard  method  of  preparing 
the  amorphous  forms  of  the  Group  V-VI  compounds  by  vacuum  evaporation  from  the 
compound  is  very  unreliable  from  the  compositional  viewpoint.  It  appears  pro- 
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bable  that  no  data  has  ever  been  reported  on  truly  Sb  Se  "  amorphous  films, 

Cm  O 

despite  considerable  literature  on  the  subject.  Several  approaches  have  been 
taken,  or  are  under  development,  which  have  yielded  amorphous  Sb2Se3  films. 

Initial  Mossbauer  studies  have  shown  that  the  Sb  atoms  in  amorphous 
Sb^Se^  (x  <  2,  y  >  3)  films  differ  in  environment  from  the  single  crystal 
form  and  therefore  appear  to  be  amorpr  ius  rather  than  microcrystalline. 


Highlights  in  the  optical  studies  to  date  are:  the  band  gap  in  Sb2Se.j 
appears  to  be  due  to  a  direct  transition;  plasma  resonance  frequencies  can 
be  accurately  determined  by  thermoreflectance  data;  photon  counter  techniques 
show  great  promise  in  determining  the  critical  points  in  band  structure? 
measurement  of  the  phase  shift  between  the  incident  radiation  and  the  photo¬ 
response  is  a  powerful  technique  for  determining  the  spectral  dependence  of  life¬ 
time  and  thus  can  yield  fundamental  data  on  differences  between  the  amorphous 
and  single  crystal  forms  of  a  material  (eg.  Se);  contacts  do  not  influence  the 
spectral  response  of  Sb^Se^  provided  that  they  are  shielded  from  the  incident 
radiation. 

We  feel  that  a  good  start  has  been  made  in  the  past  four  months  on  a 
broad  systematic  approach  toward  the  understanding  of  the  electrical  and  optical 
properties  of  the  Group  V-VI  amorphous  semiconductors.  A  better  understanding 
of  tine  basic  properties  of  these  amorphous  materials  should  lead  to  improved, 
more  reliable  and  new  device  applications. 
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THEORETICAL  RESEARCH 


The  theoretical  research  accomplished  during  the  first  four  months  of  the 
contract  period  has  been  substantial.  The  work  on  Kramers-Kronig  (Hilbert) 
transforms  is  essentially  complete  and  is  to  be  detailed  in  two  papers,  the 
first  of  which  is  now  being  submitted  for  publication  (Appendix  1).  The 
second,  which  contains  the  specifics  of  applying  the  technique  to  optical 
properties,  i.e.  to  obtaining  the  optical  functions  of  a  system  from  normal 
incidence  reflectivity  data  is  in  preparation  for  submission.  It  appears 
that  the  method  of  carrying  out  Kramers -KrBnig  transforms  which  has  been 
adopted  in  this  work  is  extremely  accurate  and  that  it  can  be  of  great  value 
for  the  analysis  of  data  as  well  as  for  obtaining  Hilbert  transforms  which 
are  of  interest  in  applied  solid-state  theory  such  as  those  of  densities  of 
states  calculated  by  Monte-Carlo  techniques.  These  are  of  interest  in  im¬ 
purity  problems. 

There  is  underway  another  effort  dealing  with  analysis  of  optical  data. 

In  this  instance  we  are  attempting  to  construct  computational  techniques  to 
extract  the  optical  constants  of  thin  films  upon  non-absorbing  substrates 
from  normal-incidence  reflectivity  and  transmissivity  measurements.  This 
is  f:  problem  which  has,  in  our  opinion,  never  been  satisfactorily  resolved 
although  much  effort  has  been  expended  toward  its  solution.  It  appears  that 
measurements  upon  at  least  two  systems  of  differing  thickness  and  a  fitting 
technique  which  includes  interference  exactly  in  some  spectral  regions  is 
necessary  for  the  proper  resolution  of  the  optical  constants  of  the  thin  film 
materials  with  which  we  are  concerned.  In  other  spectral  regions  it  appears 
possible  to  neglect  interference  and  to  sum  up  the  multiple  reflections  in 


intensity  in  order  to  find  n  and  k.  There  are,  in  addition,  spectral  regions 
in  which  interference  is  not  resolvable  and  the  exact  equations  may  be 
appropriately  averaged.  Computational  codes  for  dealing  with  this  important 
data  reduction  problem  are  in  development. 

There  is  .continuing  work  on  the  electronic  structure  of  amorphous  solids 
The  preliminary  calculation  on  Ge  indicated  that  great  care  must  be  taken  in 
the  construction  of  ensemble-averaged  structure  factors  for  amorhpous  systems 
in  order  to  eliminate  spurious  energy  eigenvalues  from  the  final  results  (i.e 
energy  eigenvalues  which  appear  in  the  solution  for  the  amorphous  system  but 
do  not  converge  to  the  proper  solutions  in  the  extreme  crystalline  limit). 
Presently  the  problem  is  being  reformulated  in  the  Wannier  (real -space) 
representation  in  terms  of  the  real-space  pair-correlation  function.  The 
same  techniques  for  averaging  computationally  over  an  ensemble  of  lattice 
configurations  are  being  used.  It  is  expected  that  this  reformulation  of  the 
problem  will  lead  to  proper  behaviour  of  the  solutions  in  the  extreme  limits 
of  the  crystalline  and  random  lattices. 


MATERIALS  PREPARATION 


A)  Single  Crystals 

No  difficulties  have  been  experienced  with  growing  large  single  crystals 
of  Sb2Te3  and  Sb2Sej.  /.  simple  horizontal  zone-refiner  with  background 
heating  to  prevent  dissociation  of  the  compound  during  melting  is  being  used 
to  grow  the  crystals.  The  shape  of  the  freezing  interface  can  be  seen  and 
controlled  for  self-seeded  single  crystal  growth. 

Growth  of  Sb2$3  has  proved  to  be  more  difficult  due  to  the  high  dissociation 
vapor  pressure  of  S  over  the  melt.  A  vertical  zone-refiner  and  a  vertical 
Bridgman  furnace  with  an  extremely  sharp  temperature  gradient  have  been 
constructed.  It  is  expected  that  one  of  these  will  prove  suitable  for  the 
growth  of  single  crystal  Sb2Sg. 

B)  Amorphous  Thin  Films 

An  amorphous  thin  film  can  be  readily  prepared  by  vacuum  evaporation  from 
a  source  of  Sb2Se3  on  t0  a  room-temperature  substrate.  This  is  the  standard 
procedure  reported  in  the  literature.  However,  we  have  found  by  chemical  and 
x-ray  emission  analysis  that  the  resulting  compositions  of  the  films  deviate 
considerably  from  stoichiometry,  generally  on  the  Se-rich  side  and  with  the 
composition  depending  greatly  on  the  temperature  of  the  evaporation  boat.  Since 
there  is  considerable  data  in  the  literature  on  so-called  amorphous  "Sb^Se^" 
films,  it  is  of  considerable  interest,  therefore,  to  examine  the  optical 
properties  of  amorphous  films  of  the  Sb-Se  system  as  a  function  of  Sb  to  Se 
ratio. 

The  type  of  apparatus  used  for  our  evaporation  experiments  was  developed 
prior  to  the  contract  and  is  generally  described  in  the  attached  publication 
(Appendix  2)  with  the  main  difference  being  that  the  substrate  temperature  can 
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be  controlled  down  to  liq.  N2  temperatures  by  insertion  of  a  copper  cold 
-finger  into  the  system. 

We  have  found  that  Se  is  preferentially  evaporated  at  low  source 
temperatures  -  particularly  below  the  melting  point  of  the  source  material 
(mp.  of  Sb2Se3  ^  610°C).  Films  approaching  stoichiometry  have  been  made  by 
raising  the  source  temperature  to  700°C,  i.e.  well  above  the  melting  point 
and  well  above  that  required  for  fast  rates  of  depositions;  but  the  films, 
although  less  Se  rich  than  at  lower  temperatures,  were  still  in  the  range 
SbxSey  (where  x  «  y).  The  results  suggest  that  Sb2Se3  does  not  evaporate 
as  a  molecular  species.  However,  the  lack  of  stoichiometry  in  the  films  may 
be  partly  a  sticking  coefficient  problem;  thus  a  comparison  of  composition  as 
a  function  of  substrate  temperature  will  be  tried  and  should  help  elucidate 
this  effect. 

Currently,  experiments  are  underway  to  use  large  evaporation  sources 
rich  in  Sb,  eg.  SbSe,  and  to  examine  the  stoichiometry  of  the  films  as  a 
function  of  subsequent  depositions  keeping  the  source  temperature  fixed  at 
^650°C.  Initial  results  show  some  success  by  this  approach.  A  thick  film 
(^20p)  of  Sb2Se3  composition  has  been  prepared  and  will  be  used  in  the  Moss- 
bauer  experiments. 

An  alternative  approach  to  obtain  stoichiometric  films  is  being  tried, 
namely,  coevaporation  of  Sb  and  Se  in  an  ion-pumped  bell -jar  system  with  the 
relative  rates  of  evaporation  being  controlled  by  quartz-crystal  oscillators. 
Early  difficulties  in  rate  fluctions  appear  now  to  be  resolved  and  again 
early  indications  are  that  films  of  controlled  stoichiometry  can  be  prepared 
by  this  technique. 

A  third  approach  to  be  tried  is  to  prepare  amorphous  films  by  sputtering, 
and  a  system  has  been  designed  to  accomplish  this.  Correspondence  between 
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source  and  deposited  film  is  expected  by  first  allowing  equilibrium  conditions 
to  be  maintained  before  exposing  the  substrate  to  deposition  (1),(2),(3). 

R.  F.  Sputtering  will  be  employed  to  clean  the  fused  quartz  substrates  before 
deposition  of  the  film. 

C)  X-ray  Emission  Analysis 

A  Norelco  Vacuum  X-ray  Spectrograph  with  Li F  analyser  and  Cu-target 
has  been  employed  to  analyse  the  composition  of  our  evaporated  films.  The 

intensities  of  the  K  ,  and  K  o  lines  of  Sb  and  Se  were  measured  to  determine 

al  ct  c 

the  Sb:Se  ratio  in  the  films.  This  is  not  a  reliable  absolute  method  of 
analysis  and  therefore  it  is  essential  to  prepare  standards  for  analysis  and, 
because  of  reabsorption  of  the  characteristic  emmitted  x-rays.  It  is  essential 
that  these  standards  are  of  a  form  closely  approximating  to  that  of  the  material 
to  be  analyzed. 

Standards  for  analysis  were  prepared  by  heating  weighed  amounts  of  Sb 
and  Se  in  a  flat-bottomed  evacuated  quartz  capsule  placed  in  a  rocking  furnace. 
The  charge  was  quenched  onto  the  flat-bottom  section  and  thus  approximated  in 
geometry  to  the  thin  evaporated  film  samples.  Despite  extreme  quenching 
conditions  (i.e.  quenching  into  a  super-saturated  salt-water  solution  held  near 
solid  CO2  temperatures)  the  resulting  discs  of  SbxSey  compositions  were  always 
crystalline  showing  that  vapor  quenching  is  much  more  effective  in  producing 
amorphous  material.  Since  the  x-ray  emission  intensities  are  thickness 
dependent,  reliable  analysis  can  only  be  obtained  on  very  thick  or  very  thin 
samples . 
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STRUCTURAL  INVESTIGATIONS 


A)  Mossbauer  Studies 

Mossbauer  studies  using  a  high  intensity  Sb  source  have  been  carried  out* 

on  two  powdered  samples  of  single-crystal  Sb2Se3  and  on  two  amorphous  films 

of  Sb  Se  ,  vacuum  deposited  onto  fused  quartz  substrates.  Very  thick  tens 
x  y 

of  microns)  samples  of  Sb^e^  were  prepared  for  these  experiments  in  order  to 

O 

obtain  the  optimum  content  of  5  mg  of  Sb  per  cnr  for  these  measurements.  Diffi¬ 
culty  experienced  in  the  thick  films  peeling  from  the  substrate  was  overcome 
by  reducing  the  size  of  the  source  aperture  thus  slowing  the  evaporation  rate. 

The  single  crystal  samples  showed  Mossbauer  spectra  which  could  be  fitted 
only  by  postulating  two  sites  for  the  Sb  atoms,  which  is  consistent  with 
crystallographic  data.  For  the  amorphous  films  a  single  Sb  site  fit  was  ob¬ 
tained  showing  that  the  very  short  range  order  differs  significantly  in  the 
amorphous  form,  i.e.  the  films  cannot  be  regarded  as  microcrystalline.  It  should 
be  noted  that  these  films  were  of  the  conventional  amorphous  variety  i.e.  non- 
stoichiometric  composition.  A  further  thick  film  sample  has  been  prepared 
which  appears  from  x-ray  emission  data  to  be  Sb2Se3  and  this  will  be  Mossbauer 
analyzed  shortly. 

A  further  Mossbauer  study,  to  be  carried  out  in  our  laboratories,  has  been 
designed  to  clarify  the  role  of  Sn  in  As2Se3Snx-glasses  on  the  composition  de¬ 
pendence  of  electrical  conductivity  and  microhardness.  A  model  in  which  a  SnSe2~ 
type  bond  is  replaced  by  a  SnSe-type  bond,  as  the  concentration  of  tin  increases, 
has  been  postulated  by  Shkol'nikov  (4)  to  explain  hardness  and  conductivity 
changes.  Borisova  et  al  (5)  conclude  in  a  recent  Mossbauer  study  that  the  only 
Sn4+  occurs  and  that  Sn  is  bonded  only  to  Se;  however,  an  examination  of  their 
data  leads  us  to  conclude  that  a  distribution  of  ionic  states  is  present  which 


*  In  collaboration  with  Dr.  Stan  Ruby  of  Argonne  National  Laboratory. 
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could,  in  fact,  be  that  predicted  by  Shkov'nikov.  Further,  Borisova  et  al 
do  not  find  the  correlation  of  intensity  with  coordination  number  previously 
found  in  glasses. (6)  The  Mossbauer  experiments  on  ASgSe^Sn^-gl asses  will  be 
repeated  with  narrower  source  lines  and  better  statistics  to  determine  the 
degree  of  local  order  and  the  distribution  of  ionic  states  of  Sn.  A  test 
of  the  relationship  of  Mossbauer  intensity  to  coordination  will  also  be  made. 

B)  X-ray  Studies 

The  experimental  equipment  consisting  of  a  Laue  camera,  a  Weissenberg 
and  precession  camera,  in  connection  with  a  Norelco  generator,  has  been 
set  up  for  investigating  the  Sb?Te„  Se.  system.  In  addition  to  verifying 
and  possibly  refining  the  exact  structure  of  Sb,Je3  x_ray  data,  it  is  of 
particular  interest,  to  study  the  changes  in  the  structure  due  to  incorporation 
of  Se  atoms.  The  question  of  the  specific  ordering  is  of  great  importance 
since  it  has  a  strong  influence  on  the  transport  properties. 

According  to  Teramoto  and  Takayangi  (7)  in  the  Sb?Te~  Se  solid 
solutions  2/3  of  all  Te  atoms  can  be  replaced  by  Se  without  changes  of  the 
rhombohedral  (hexagonal)  structure.  Due  to  the  smaller  atomic  radius  of  Se, 
the  lattice  dimensions  of  a  and  c  decrease  with  increasing  amounts  of  Se. 

When  2/3  of  the  Te  atoms  are  replaced  by  Se,  the  lattice  constant  change 
amounts  to  about  4%  for  a  and  3 %  for  c.  The  antimony  may  go  into  either  the  c 
or  a  positions  or  be  statistically  distributed  into  both.  Since  the  atomic 
scattering  factors  of  Te  and  Se  vary  considerably,  the  scattering  amplitude 
will  differ  considerably  for  these  various  possibilities.  Once  the  correct 
model  is  established,  the  parameters  can  be  refined  by  a  least  squares  method 
using  the  intensities  of  all  reflections.  Changes  of  inter-atomic  distances 
will  be  Correlated  with  transport  properties. 
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OPTICAL  PROPERTIES 

Work  related  to  this  program  on  the  properties  and  growth  of  Group 
V-VI  compounds  was  already  underway  before  the  commencement  of  this  contract. 

In  particular,  a  large  amount  of  work  on  the  growth  of  single  crystals  of 
Sb2Te3  -  Sb2Se3  solid  solutions  together  with  the  measurements  of  their  optical 
properties,  was  nearing  completion.  However,  analysis  of  the  data  was 
carried  out  under  the  contract  and  the  results  of  this  work  have  been  accepted 
for  publication  (Appendix  3). 

A)  Reflectivity 

The  optical  reflectance  and  transmission  of  single  crystal  specimens  of 
solid  solutions  of  the  Sb2Te3  ^Se^  compounds  were  investigated  in  the  photon 
energy  range  0.04  eV  to  6.5  eV.  A  wide  range  of  solid  solutions  from  x  =  0 
to  2  were  prepared  all  having  rhombohedral  structure  and  were  found  to  exhibit 
properties  closely  corresponding  to  the  isostructural  compounds  Be2Te3,  B12Se3. 
Small  indirect  band  gaps  approaching  semimetallic  values  were  obtained  for  all 
compositions.  The  high  free-carrier  concentrations  gave  rise  to  plasma 
effects  in  the  near  infra-red.  The  chemical  bonding  concepts  applied  to  similar 
systems  seemed  applicable  with  only  slight  modification. 

A  result  of  some  significance  to  this  program  was  that  both  the  reflectivity 
and  transmission  data  show  that  the  absorption  edge  of  orthorhombic  Sb2Se3 
corresponds  to  direct  transitions,  at  least  for  non-polarlzed  light  perpendicular 
to  the  cleavage  plane,  i.e.  with  the  electric  vector  parallel  to  the  (a)  and  (c) 
directions.  This  is  In  contrast  to  the  assignment  of  an  indirect  transition  by 
Shutov  et.  al .  (8)  and  is  of  utmost  importance  for  future  comparisons  of  the 
band  structure  of  Sb2Se3  In  the  amorphous  and  crystalline  form. 

Experiments  are  currently  underway  to  measure  the  reflectivity  and 
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transmission  data  as  a  function  of  polarized  light  In  the  three  primary 
crystallographic  orientations  of  Sb2Se.j.  Some  difficulty  has  been  experienced 
in  preparing  good  reflectivity  surfaces  cut  perpendicular  to  the  cleavage 
plane,  or  thin  optical  transmission  samples  of  the  same  orientation,  due 
to  the  propensity  of  the  crystals  to  cleave  easily.  An  encapsulation  tech¬ 
nique  In  a  polyester  casting  resin  shows  some  promise  of  overcoming  th4s 
difficulty.  The  suitability  of  various  chemical  etches  for  preparing  surfaces 
for  reflectivity  measurements  has  also  been  investigated.  For  Sb2Se.j..  2 
parts  HC1 :  1  part  H202  by  volume  has  been  found  to  produce  a  satisfactory  surface; 
and  for  Sb2Te3»  the  best  etch  appears  to  consist  of  14  parts  Acetic  Acid 
(saturated  with  Citric  Acid)  1  part  Bromine. 

Although  no  difficulties  have  been  experienced  In  preparing  good  amorphous 
thin  film  specimens  for  optical  measurement,  work  has  been  held  up  In  this  area 
mainly  by  lack  of  control  of  stoichiometry  of  the  films  (see  Preparative  Section) 
and  by  lack  of  a  suitable  analytical  techniques  to  accurately  obtain  the  n  and  k 
values  for  comparison  with  data  on  single  crystals  (see  Theory  Section^. 

B)  Thermally  modulated  reflectivity 

The  thermally  modulated  reflectivity  of  single  crystalline  materials 
Sb„Te,  Se  (x  =  0  -  2)  was  measured  in  the  UV,  visible  and  IR  region. 

Experiments  were  performed  In  the  IR  region  (8-28m)  in  order  to  confirm  the 
existence  of  plasma  resonance  of  free  holes.  The  reflectivity  structure  due 
to  this  resonance  was  rather  broad  due  to  strong  free-carrlers  absorption  and 
Interband  (light  and  heavy  hole  bands)  transitions,  which  mask  the  resonance 
effect.  However,  the  thermal  modulation  gave  precisely  the  resonance  fre¬ 
quencies  Up.  The  knowledge  of  this  resonance  frequency  Up  together  with  the 
Hall  coefficient  and  resistivity  can  determine  the  optical  effective  mass  and 
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the  relaxation  time  due  to  free-carrier-phonon  scattering.  This  data  has  been 
submitted  for  publication  (see  Appendix  4). 

The  thermomodulated  reflectivity  was  also  performed  for  SbgTe^  in  the 
visible  and  near  UV  region  using  a  modified  Cary  14R  spectrophotometer  to 
confirm  the  existence  of  critical  points  corresponding  to  transitions  at  points 
of  high  symmetry  in  the  Brillouin  zone.  Preliminary  results  yield  quite  sharp 
structure  which  indicates  agreement  with  identification  of  these  critical  points 
in  the  work  of  Sobolev  et  al .  (9)  from  low  temperature  reflectivity  measurements. 
The  structure  of  the  thermomodulation  spectrum  agrees  with  the  photon  counter 
(see  Transmission  measurements). 

So  far  it  has  been  found  difficult  to  extend  measurements  beyond  the 

region  3000A  to  7000A°  due  to  cut-off  of  the  photomultiplier  and  to  low  light 

intensity  in  the  UV  region.  These  problems  are  expected  to  be  solved  using 

the  RCA  C31015C  "Quantacon"  photomultiplier  having  an  extended  response  from 

2000  to  9000A°  and  using  the  McPherson  225  Vacuum  UV  monochromator.  The 

Sb  Te  Se  compounds  are  very  suitable  for  the  thermomodulated  spectroscopy 
l  0-x  x 

for  x  =  0  to  2  due  to  their  low  resistivity.  However,  for  x  >  2,  the  thermo¬ 
modulation  is  difficult  as  the  resistivity  increases  with  x  and  other  forms 
of  modulations  have  to  be  applied.  Electric  field  and/or  photomodulation  seem 
to  be  the  most  appropriate. 

C)  Transmission  measurements 

Direct  transmission  measurements  were  performed  on  thinly  cleaved  specimens 

of  Sb  Te  Se  (x  *  0  to  3)  in  the  2000  to  10,000A°  wavelength  region.  This  • 

2  3-x  x 

photon  energy  range  is  higher  than  the  direct  gap,  and  because  the  absorption 
coefficient  is  high  (^10^  to  10^  cm-*)  in  this  range,  a  sensitive  detector  must 
be  used  even  for  very  thin  samples  (^lu)-  Also  stringent  precautions  must  be 
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taken  to  avoid  any  stray  light.  The  latter  problem  was  successfully  solved 
using  an  0-rlng  sealed  vacuum-tight  blackened  tube  with  the  sample  sealed 
Inside.  Using  a  photon  counter,  a  transmission  of  10"^  to  10"^  was  measured 
with  a  S/N  ratio  of  100/1  to  100/3  provided  good  stability  of  the  light 
source  was  maintained.  This  set-up  Is  expected  to  be  Improved  using  a  double 
beam  method,  l.e.  a  double-channel  photon-counter  with  an  on-line  computer. 

The  transmission  spectrum  showed  fine  structure  corresponding  to  critical 
points  In  the  Srlllouin  zone  In  agreement  with  published  data  for  x  =  0  by  V.  V. 
Sobolev  et  al.,  (9)  and  for  x  =  3  by  S.  D.  Shutov  et  al .  (8)  Our  structure  was 
much  sharper  than  that  obtained  by  the  above  Investigators,  even  In  comparison 
with  their  low  temperature  reflectivity  data. 

A  sharp  transmission  window  at  ^^OOA0  was  shown  by  all  samples,  which 
coincides  with  sharp  fall  of  the  reflectivity.  This  effect  Is  probably  due 
to  sharp  fall -off  of  the  density  of  states,  In  agreement  with  the  band  model 
of  Isostfuctural  Sb2Te3  by  Borghese  and  Donato  (10)  rather  than  to  a  plasma 
resonance  of  electrons  In  the  valence  band.  The  latter  effect  Is  expected  to 
occur  at  a  higher  energy  range  (^10  -  14  eV)  and  can  provide  the  coordination 
number  and  the  degree  of  loniclty  of  the  bonds. 

D)  Photoemission 

Photoemission  measurements  can  yield  the  electronic  density  of  states, 
the  electron-electron  scattering  length  and  the  conservation  of  the  electron 
momentum  can  be  obtained.  The  latter  Is  of  considerable  Importance  for  the 
study  of  the  difference  between  amorphous  and  single  crystalline  materials  as 
far  as  the  selection  rules  and  the  degree  of  localization  of  electronic  states 
are  concerned. 

A  special  phototube  has  been  designed  and  Is  under  construction  which 
will  allow  measurements  on  freshly  cleaved  surfaces  of  single  crystals  or  on 


amorphous  films  evaporated  In  the  ultra-high  vacuum  of  10“®  to  lO"*1  torr. 

A  vacuum  UV  Monochromator  (McPherson  225)  will  be  used  as  the  radiation  source. 

This  monochromator  has  been  modified  by  attaching  a  Varlan  4001/sec  noble  gas 
pump  to  assure  oil-free  operation,  which  Is  a  necessary  condition  for  the 
windowless  operation  of  the  light  source  for  energies  over  12  eV.  The 
possibility  of  utilizing  the  radiation  from  the  Synchrotron  at  the  University 
of  Wisconsin  Is  also  being  explored. 

An  electronic  circuit  has  been  constructed  to  measure  the  photoemission 
current  10"11  AtolO-^A).  The  main  part  consists  of  a  feedback  Kelthley 
604  or  640  electrometer,  and  a  Kelthley  82  and  a  PAR  HR  -  8  Lock-In  amplifier 
together  with  a  bridge  circuit  for  rejection  of  the  capacitive  component  of 
the  photoemission  current. 

E)  Photoconductivity 

Photoconductivity  measurements  have  been  performed  with  emphasis  on  optical 
and  transport  properties.  An  experimental  technique  of  measuring  the  spectral 
dependence  of  the  lifetime  t  of  photogenerated  carriers  was  developed  which 
employs  a  measurement  of  the  phase  shift  >p  between  the  square  light  pulses  and 
the  first  harmolc  frequency  of  the  photocurrent  response  using  a  Lock-in  amplifier. 

This  method  was  applied  first  to  the  rather  well -Investigated  element  Se  in 
the  amorphous  and  single  crystal  (trigonal)  form  at  90°K  and  300°K.  It  was 
found  that  the  photoconductivity  of  single  crystal  of  trigonal  Se  comprises  both 
valence  to  Conduction  band  transitions  as  well  as  transitions  from  states  inside 
the  forbidden  gap  Into  the  conduction  band.  In  the  amorphous  Se,  only  the 
valence  band  to  conduction  band  transitions  contribute  to  the  photoconductivity. 

In  fact,  in  amorphous  Se  the  photoconductivity  onset  occurs  at  higher  energies 
than  the  bend  gap,  which  has  been  explained  (11)  by  strong  exclton  absorption. 

This  work  has  been  submitted  for  publication  (Appendix  5). 
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Similar  results  were  obtained  for  amorphous  films  of  Sb£Se3  with  an 

excess  of  Se.  While  the  photoconductivity  onset  for  these  films  occured  either 

at  0.8y  or  at  l.Op,  the  absorption  edge  Is  located  at  longer  wavelengths  (  1.2 

to  1.3m)*  Further  experiments  on  the  photoconductivity  of  stoichiometric  amor- 

phous  films  of  Sb  Se  should  Indicate  whether  the  above  effect  Is  due  to  exclton 
2  3 

absorption,  similar  to  case  of  amorphous  Se  films,  or  the  light  scattering  In 
amorphous  films  or  to  optical  absorption  due  to  excess  Se  atoms  In  non-stolchlo- 
metrlc  films. 

The  photoconductivity  of  Sb^Se^  films  was  measured  using  evaporated  Au  or 
Sb  contracts.  The  overall  spectral  response  was  the  same,  for  both  contacts  con¬ 
trary  to  the  work  of  Lyubln  et  al .  (12)  who  maintain  that  Au  contacts  were  block¬ 
ing  and  Sb  Injecting,  and  that  they  produce  a  different  photoconductlve  spectral 
response.  Our  optical  transmission  measurements  on  Sb2Se3  films  show  an  absorp¬ 
tion  peak  at  0.75m  In  accordance  with  (13)  and  sometimes  additional  structure  at 
1.0m  and  0.85m.  Results  In  (12)  may  be  due  to  the  use  of  a  sandwich  structure 
with  unavoidable  Illumination  of  contacts  producing  photovoltaic  effects,  while 
In  our  case  contacts  were  carefully  shielded  from  Illumination. 

The  photoconductivity  In  orthorhombic  single  crystals  of  Sb  Se  was  Investi¬ 
gated  as  a  function  of  polarization  of  the  eloctrlcal  vector  E0  of  the  Incident 
light  with  respect  to  crystallographic  axes  of  the  single  crystal.  Edr.ft  was 
||  c  with  evaporated  gold  contacts.  For  this  orientation,  the  photocurrent 
showed  a  maximum  at  1 .05y  and  a  secondary  maximum  at  0.75m.  For  E  | |  a  the 
first  maximum  was  shifted  towards  1.075m  (In  accordance  with  the  shift  of  the 
absorption  edge)  while  the  secondary  maximum  remained  at  0.75m.  Also,  another 
maximum  was  found  at  0.575m.  The  secondary  maximum  at  0.75m  was  not  reported 
In  earlier  work  on  single  crystals  (13)  but  was  reported  to  occur  In  Sb2Se3 
amorphous  films  (13),  In  agreement  with  our  results. 
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TRANSPORT  PROPERTIES 

Hall  Effect  Measurements 

There  Is  some  evidence  to  Indicate  that  Inorganic  materials  with  a  chain- 
like  crystal  structure,  e.g.  single  crystals  of  selenium  (14),  exhibit  a 
barrier-type  photoconductivity.  The  phctoconductlve  mechanism  In  some  compounds 
containing  selenium,  e.g.  PbSe,  and  related  chalcogenldes,  has  long  been  In 
contention  with  the  barrier  theory  receiving  strong  consideration  (15).  Since 
the  arsenic  and  antimony  chalcogenldes  are  related  structurally  and  composition- 
ally  to  the  above  materials  It  Is  essential  that  methods  developed  tc  Investigate 
transport  and  photoconductlve  mechanisms  will  allow  one  to  distinguish  between 
carrier  generation  and  barrier  layer  theories.  A  Hall-effect  or  photo-Hall 
effect  measurement  In  which  the  frequency  of  the  applied  electric  field  can  be 
varied  up  to  high  values  (^  MHz)  will,  by  virtue  of  shorting  potential  barriers, 
allow  this  distinction  to  be  made. 

Initially,  the  Hall  coefficient  for  single  crystals  of  selenium  will  be 
measured  as  a  function  of  frequency  of  the  applied  electric  field.  The  fre¬ 
quencies  used  will  range  of  20  Hz  to  1  MHz.  The  resistivity  will  also  be 
measured  for  the  same  frequency  range. 

A  typical  6  contact  Hall  sample  will  be  used  together  with  3  RF  switches 
to  measure  current.  Hall  voltage  and  voltage  drop  due  to  resistance.  The  Input 
signal  will  be  obtained  from  a  Hewlett-Packard  oscillator  through  a  RF  trans¬ 
former  which  has  a  flat  response  from  20  Hz  to  5  MHz.  The  actual  voltage 
measurements  will  be  made  by  a  Kelthley  lock-in  amplifier  system.  Preliminary 
measurements  Indicate  no  distortion  of  wave  form  up  to  1  MHz  and  signals  as  small 
as  10y  vOlts  have  been  measured  to  an  accuracy  of  a  few  percent. 
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ABSTRACT 


It  Is  snown  that  the  Hermite  functions  form  a  compact,  swiftly 
convergent  representation  of  the  Hilbert  transform  operator  (Kramers- 
Kronig  transform).  Application  is  made  to  the  practically  Important 
case  of  the  Hilbert  transform  of  a  function  non-zero  only  over  a  finite 
region.  Hermite  functions  form  a  natural  expansion  set  in  such  cases. 

If  the  function  to  be  Hilbert  transformed  has  discontinuous  derivatives 
(critical  point  structure)  at  a  set  of  values  u^,  then  the  transformed 
function  also  has  discontinuous  derivatives  at  the  same  set  of  w^.  The 
number  and  u>  position  of  critical  point  structure  is  shown  to  be  invariant 
to  Hilbert  transformation.  The  practical  significance  of  measuring 
experimental  optical  data  only  at  the  zeros  of  Hermite  polynomials  is 
stressed.  The  zeros  and  weight  functions  for  Hermite  integration  for 
n  =  20,  26,  50,  76,  150,  and  300  is  given  in  an  appendix. 


I.  INTRODUCTION 


Dispersion  relations,  starting  with  the  classical  work  of 
1  2 

Kronig  and  Kramers,  have  been  applied  to  an  ever  widening 

VS 

variety  of  problems  in  physics.  Elegant  formalisms^  now 
exist  so  that  the  mathematical  structure  of  the  Hilbert  transform 
Is  fully,  and  completely  understood.  The  practical  application 
of  the  Hilbert  transform  to  physics,  especially  to  the  analysis 
of  experimental  data,  has  been  problematic  for  two  principle  reasons 
(1)  Direct  integration  of  (necessarily)  discontinuous  data  with 
singular  kernels  is  computationally  impossible;  (2)  Data  do  not 
exist  over  the  whole  real  line  (-  ®  <  ») ,  but  only  over  a 

finite  range  or,  at  best,  a  sequence  of  overlapping  finite  ranges. 

The  first  problem  has  generally  been  treated  by  "smoothing"  the 
data  over  a  region.  If  the  smoothing  was  achieved  by  expanding 
M  data  points  In  a  Taylor  series  to  order  N  (M  >  N) ,  then  the 
smooth  function  Is  continuous  In  its  zeroth  through  Nth  deriva¬ 
tives  In  the  region,  but  discontinuous  in  the  sense  that  N  +  1th 
and  higher  derivatives  are  set  equal  to  zero  (Removing  informa¬ 
tion  from  the  Hilbert  transform. )  At  the  region  boundaries  all 
derivatives  are  discontinuous.  The  effect  of  these  discontinuities 
and  forced  zeros  on  the  resultant  Hilbert  transform  Is  generally 
unknown.  This  problem  is  usually  solved  by  making  plausibility 
arguments  that  if'  the  transform  function  is  smooth  and  or  If 
applying  the  Hilbert  transform  twice  returns  approximately  the 
negative  of  the  original  smoothed  function,  then  the  transformation 
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was  successful. 

The  second  problem  exacerbates  the  first.  Generally  the 

data  are  continued  above  the  last  data  point  by  means  of  a 

Laurent  expansion  ( ~  £  a^  u^11  ^n),  where  an  are  expansion 

n 

coefficients,  and  ^  is  the  measured  independent  variable.  Below 
the  first  data  point  a  Taylors'  series  is  the  expansion  set.  Again 
plausibility  arguments  are  used.  (We  neglect  completely  in 
this  paper  the  problem  of  extrapolation  in  cases  where  significant 
structure  exists  in  an  unmeasured  region. ) 

The  increased  application  of  Monte  Carlo  methods  to  physical 
models  yields  numerical  results  which  cannot  be  represented  in 
closed  form.  Thus  the  calculation  of  the  Hilbert  transform  of 
such  results  is  also  difficult  as,  for  example,  in  the  important 
case  of  finding  the  Hilbert  transform  of  a  density  of  states. 

In  this  paper  we  consider  a  mathematically  exact  treatment 
of  this  problem  which  is  simple  to  use,  can  be  applied  to  a 
variety  of  cases  ,  automatically  least-square  fits  discontinuous 
data  over  the  entire  measured  region,  provides  a  natural  form 
of  extrapolation  function,  and  in  the  special  case  of  optical 
data  allows  a  simultaneous  best  fit  to  all  of  the  variables 
fig  (  uj)  >  e2(u))>  InR  (u>)>  cp  ( u>)  N(  u>)  ,  k(  u>) ,  and  T(u>)  in  terms  of 
one  expansion  set,  where  and  e2  an^  N(  u>)  and  k(  «>)  are  thts.  real 
and  imaginary  parts  of  the  dielectric  function  and  the  index 
of  refraction  respectively,  R  is  the  normal  ref lectlvLty ,  the 
phase,  and  T  the  fractional  transmitted  intensity.  In  Section 
II  we  discuss  the  formal  derivation  of  our  technique.  In 
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Section  III  we  give  the  results  of  two  applications:  the  Hilbert  transform  of  the 
density  of  states  of  Monte  Carlo  calculations  of  a  simple  cubic  s  band  and  face 
centered  cubic  palladium.  Section  IV  summarizes  our  results  and  discusses 
application  to  other  systems.  The  numerical  application  to  the  analysis  of  normal 
reflectivity  data  will  be  detailed  In  a  further  publication. 


II.  FORMALISM 


A.  The  Hilbert  Operator 


As  is  well  known,  the  Kramers-Kronig  relation  between  the 
real  and  Imaginary  parts  of  a  response  function  arise  because  of 
causality.  If,  in  the  long  wave  length  limit,  a  time  varying 
field  K(t)  Is  Impressed  upon  a  media  at  t  «  0,  then  the  linear 
response  f?(  t)  must  also  be  zero  for  t  <  (),  and  is  given  by 


D(t)  e(t)  E(T)  (1) 

where  e  is  the  generalized  susceptibility.  Upon  taking  the 
complex  Fourier  transform  of  (1)  we  obtain  the  Kramers-Kronig 
relations : 


|>  (.  cs(u)'  )<lc' 
*1<">  -  1  n  .1  — : - 


to 1  -  u> 


,  .  p  p00  (  €  i  (  ul)  “  l)*l 

Co(  (o)  =  ”  —  - 

n  ... '  - 


to' 


U)  -  0) 
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»h*re  «,  and  *  denote  the  real  and  Imaginary  parte  and  the  1 
comes  from  the  fact  that  =  1  Tf 

,  ,  l[  J  1‘  If  we  reSard  the  principle 

relation  ^e^ralS  ^  ^  ?  «"  Motional 


(‘1  -  1)  =  P  («2) 
*2  "  "  P  (ej  -  1)  ‘ 


(3) 


or  that 


(«i  - 1)  -  -  p2  (Cl .  i) 


o) 

»Meh  makes  manifest  that  P  i.  an  antlunltary  operator. 

pair  "rrticians  the  reiawons  (2)  f°™ a — — 

P  ir.  The  cloee  relationship  between  Hilbert  and  Fourier 
transformations  is  seen  in  that-  t-h«  * 

tlon  1  ,  hat  the  kernel  for  Fourier  transform*- 

tion  is  given  by 


kp  ("  ■  •*)  r  **  «<>•  (.  -  w-jt 

o 

"hereae  the  kernel  for  the  Hilbert  transformed 


(5) 


on  is 


%  (.  -  *')  «  *£dt  sin  („  -  ,;)t  . 
A  more  familiar  form  for  (5)  and  (6)  is 

KF  («  "  w1)  -  fi  ( a)  -  a, » ) 


(6) 


(7) 
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(a  Dirac  delta  function)  and 

/ 

Kjj  (a)  -  <!)')  =  TT  - — - >  (8) 

(  0)  "  OJ '  ) 

which  we  define  as 

U  -  »')  s  P  .(»  -  »')  (9) 

The  relations  (5)  and  (6)  form  representations  of  the 

operators  K„  and  Ku.  Any  denumerable,  complete  set  of  functions 
r  n 

defined  over  the  whole  real  line  may  be  used  to  form  a  represen¬ 
tation  for  (5)  since  in  terms  of  any  complete  orthonormal  set 
*n(w)  we  have, 

6  ( <»>  -  co '  )  -  E  \jrn  ( to)  \|rn  ( to'  )  ,  (10) 

whore  the  sum  extends  over  all  of  the  members  of  the  set. 

An  analogous  representation  for  the  operator  P  ( m  -  co ' )  may  be 
constructed  as  follows: 

From  (6),  or  by  direct  integration  we  have  that  p  is 
given  by 

p  («-«')“  ^  r * c,t  sgn  (fc) el(u> " u,,)t  (ii) 

TT  J 

—  oo 

where  CON  ( t)  is  the  signature  or  signum  function  (-1  for  t  <  0 
and  +1  for  t  >  0  and  0,  say,  if  t  -  0) .  We  recall  that  the 
3.1gnum  function,  like  the  delta  function,  Is  a  generalized 
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MMM 


function 


4 


and  has  the  property  that 


SON  (x  -  x1)  =  2  s  (x  -  X') 


(12) 


Equation  (11)  allows  us  to  construct  an  infinite  number  of 
representation  of  the  Hilbert  operator  by  simply  finding  the 
representation  for  exp  (  i<ut  )  in  terms  of  the  complete  set  and 
using  (11).  Here,  and  for  the  rest  of  this  paper,  we  will  con¬ 
centrate  on  a  single  representation,  the  Hermite  functions, 
tn(o))*  We  recall'  that  the  nth  order  orthonormal  Hermite  function 
is  found  from  the  corresponding  Hermite  polynomial  Hn(aO  by 


2 

-* 


»„(•>) 


(13) 


where  Nn  is  a  normalization  factor  given  by 


=  (2n  n!  V^)'1  •  (14) 

The  most  important  property  of  the  Hermite  functions  we  need  is 
that  up  to  a  trivial  factor  in,  the  Fourier  transform  of  a 
Hermite  function  of  order  n,  is  the  same  Hermite  function: 


(15) 


Multiply  both  sides  of  Eq.  (15)  by  ^f)(  m) , 
we  have  a  representation  for  e*u,S 


and  utilizing  eq.  (10), 
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eiu,t  =  -\j2n  E  ln  ^n(t)  ^n(  w)  ; 
n 


and  o 


(16) 


s;  (-1)™  ,m(t)  tmU<) 


m 


( 17) 


placing  (16)  and  (17)  In  (11)  we  have  that 


P(u>  -  a)')  =  (-1)  E  E  (I)”  (-l)m  irn(  u>)  (u)')  S'm 


n  m 


(18) 


where  S'  la  given  by 


Sn,n“  1“  dt  S™  (*)  ♦„  <*>  ' 


(18) 


We  notice  that  S'  forms  the  representation  of  the  s.tgnum 
operator  In  terms  of  the  Hermite  functions.  We  recall  that  the 
Hermlte  functions  are  simply  odd  or  even  by  the  rule:  ln(-t)  = 
(-l)n  ^n(t)*  The  Integral  for  negative  t  may  be  reduced  to 


Snm  =  ^“1)m  +  "  +  1  I  dp 


(19) 


50  that  the  whole  integral  is  given  by; 


SU~  t1  +  (-D 


m  f  n  +  1 


1  I  dl  k  o  Ue 


(20) 


'c.  (20)  makes  clear  that  S'  ia  aero  unless  n  h  m  is  odd.  Let 
n  assume  that  this  is  true.  Then  (20)  may  lie  written 
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Rquatlon  (18)  makes  clear  that  s  is  a  symmetric,  real  matrix 
so  wo  may  assume  m  >  n  .  We  recall  that  the  gamma  function 
(r(ni-l)  «-•  nl)  has  the  property 


i’(2%)  -,L.2p,/-l/2  r(y,)  r( v,  v  1/2)  . 

\  8rr 

The  ratio  of  the  factorials  involving  n  -j-  m  reduce  to 


(26) 


r(mfn-Pr)  1  0(m»-n-2r)-l/2  „  ,  riH-n-2rN 

r(STH^-  2  r  <-tH 


(27) 


Let  m  -  n-i-2h-|-l,  then  after  simplifying,  we  find 


m  fc  n.1  (n  |~2h  |-1)1  ,  . 

’n,  ni-2h|  l  ~  '/  n  '  ' 


x  E 
r=  0 


(-  !.)  ( 2n<-2h-2r-l)  J ! 
r !  ( n+2h- r+1 ) ! ( n-r) ! 


(28) 


The  latter  relation,  although  correct  In  principle  to  all  orders,  becomes 
Inaccurate  for  purposes  of  computlon  above  n,m  >35.  Therefore,  Snm  was 
found  from  (1ft)  by  simply  Integrating  using  the  zeros  of  the  500tn  order 
Hermlte  functions  as  outlined  In  Appendix  A.  The  latter  technique  Is 
extremely  fast  and  up  to  order  n+m  -•  500,  Is  exact 


Returning  to  our  original  relation  (1.8),  wo  havG 


0)  -  0) 


')  =  (-1)  E  E  (J)”  (-i)’"'Sh|-1  ,  ,  ,  ,v 

n  h  viv  wf  *n-i-?h+l  '  ■  ) 


n,  m-2h  i-l 


(29) 


Let  us  define  S  such  that 
n,m 


l’(oj  -  a)' )  =  '»  .  •;•.  ( oi)  •)  ( «,»,) 

nm 


(30) 
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Thcn  the  relation  between  S  and  S'  with  m  -  nn-2h+l  is  given  by 


Snm  =  ;  Smn  “  '  Snm 


(31) 


The  compactness  of  the  operator  P  can  be  seen  either  from 
equation  (4)  or  from  the  fact  that  for  “fixed  n 


lim  S  „l0,  , ,  -♦  0 

v,  .  n  f  n-(-2h+l 


(32) 


2h  —  i 

at  least  as  fast  as (2  ( 2h  +  1) 1 1 )“  .  The  latter  is  easily  proved  by 

direct  expansion  of  equation  (28). 

Wo  now  take  up  the  problem  of  Hilbert  transforms  of  functions. 
This  problem  is  especially  of  in. -crest  in  exploring  the  Inter¬ 
relationship  of  optical  proper-ties. 


B.  Optical  properties 


We  have  seen  in  the  last  subsection  that  the  Hermite  functions 
form  a  compact  representation  of  the  Hilbert  transform  operator. 
Here  we  will  exploit  the  fact  that  those  functions  also  possess 

simple  properties  In  the  complex  plane.  Two  such  fund; ions  are 

1 

of  practical  interest:  The  complex  square  root  and  the 

homographic  transformation.  The  first  gives  the  relationship 
between  the  dielectric  function  and  the  index  of  refraction; 

The  second  the  relation  between  the.  log  of  the  reflectivity  and 
the  index  of  refraction.  oUr  purpose  hero  is  ‘o  show  that  a 
single  set  of  expansion  coefficients  of  the  II'-. :r.i '  fund. ions 
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eharacterizcs  all  the  optical  functions.  In  particular  we  will 
show  the  inter-  re  la  tions  hips  between  the  expansion  coefficients. 
Let  us  assume  that  we  have 


and 

N  =  n  +  ik  =  s  U)  (34) 

n 

where  C  ,  dm  are  complex  expansion  coefficients  which  factor 
into  real  and  wholly  imaginary  parts  as  n,m  increase  (decrease) 
by  1.  Then  we  have  N  =  (e)1^2  or  that 


Cm  E  dnd£  J  cHn((o)  *n  (to)  ♦/») 

n  i  ,  —co 


_  ,  ,  Tm 

~  n i  n  *  n* 


(35) 


Again  exploiting  the  extremely  simple  "Clebsh-Gordon"  coefficients 
for  the  Hex  •ml  to  functions,  oq.  (23),  we  have  after  a  little 
algebra,  and  assuming  that  i  >  n 


“  Nn  N»>  Via  2P  O  <£> 


i 


H 


n-i-i-2r 


<V§ "» <v§ x)  ,rxa  ax 


(36) 


I 
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Kxpanding  the  latter  again  using  eq.  (23)  we  have 


N 


n 


x 


_L_ 

V* 


N, 


m 


2 


r 


rl 


min(  m,n+4-2r)  . 

E  2n  hi 

h=0 


(37) 


/mwn+  i-2rw  lx 
'h'^  h 


n+nW-i-2r-2h 


h 


x  Hmnvi-Ji-2r-21i 


(0) 


which  makes  clear  that  m  has  the  same  parity  symmetry  a3 
n  i-  Let  us  assume  that  n  +  m  +  i  is  even.  Then  we  have 


Nm  Nn 


min(m,n-i-/-2r) 

£ 

h=0 


2rhh  rl 


hi 


n  t  in  nl  .  oT,  i  mmi-£-2r-2h 

<")  (f.)  <;•’)  <  "hr  >  (^)  — s — 


nirni  £-2r-2h  ,  o  ™  \ 

f  _1  \ - 75 - (  m-nH-;/-2r-2h)  1 

'  1  2  /  • 


(38) 


which  shows  that  l"’  Is  also  a 

I J  £ 

the  utilization  oi’  eq.  ( 32)  i», 
stantially  easier  *  nan  t  he  I’orm 
suggest.  We  will  i-ium  •  a 

Finally  we  con.-.M-r  \k-\U-  . 


compact  operator.  We  stress  that 
practical  applications  is  sub- 
oV  eq.  (38)  might  initially 
ther/.e  again  In  gee;  ion  ITT. 


‘  : i« •  homograph  1  e  •  ra 


i’<  M-ma  t.  ]  on 
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(39) 


and  the  related  trams  format  Jon 


in  %  =  |n( a'£  -  p)  -  |n (yz  +  6)  (40) 

whore  here  a  -  3=7=  6  =  1,  Z  is  the  complex  reflectivity 
Reicp  and  z  the  complex  index  of  refraction.  Since  |Z  I ,  the 
reflectivity,  is  bounded  by  1,  we  have  that  jgn(Z)  =  fln(  1  +  (Z-l)) 
and 


Exploiting  equation  (4l)  on  both  sides  of  eq.  (40),  and  by 
expanding  both  Z  and  a. Into  Hermite  functions,  and  making  repeated 
use  of  oq.  (23),  wo  would  come  to  the  conclusion  that  the  operator 
connecting  the  two  sides  of  eq.  (40)  would  also  be  compact.  How¬ 
ever  we  believe  the  analog  of  equation  (38)  for  the  logarithm  is 
too  complicated  for  practical  application.  We  conclude 


that  if  a  sufficiently  great  number  N  were  chosen,  then  an 
expansion  of,  for  example ,  any  one  of  the  optical  functions  in 
the  first  N  Hermite  functions  with  <  •xpam;  I  on  coefficients  C 
guarantees  that  all  of  the  roiatf  i  fun<m  Ions  may  be  found  with 
an  error  no  larger  than  any  rlv-n  small,  nun-bop  e. 
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III.  RESULTS  AND  APPLICATIONS 


Our  discussion  in  the  last  section  was  purely  formal  and 
essentially  hinged  on  showing  the  compactness  of  various  operators, 
in  terms  of  Hermite  functions.  Here  we  will  be  entirely  concerned 
with  the  practical  application  of  our  formal  results  to  real  sys¬ 
tems.  We  will  also  slightly  extend  our  definition  of  expansion 
and  furthermore  show  compactness  in  the  important  case  of  finite 
data  ranges. 


A.  s-band  Density  of  States 

First  we  consider  a  problem  in  which  the  function  to  be 
transformed  is  defined  everywhere.  We  apply  our  results  to  the 
Important  problem  of  deriving  the  "F"  function  for  the  density 
of  states  of  s  states  in  the  simple  cubic  lattice: 


F(  oj)  =  ^  f  "  du,' 

tt  J  0)  -  U) 


(42) 


where 


g(  co * )  =  —^5  f  dh6  (E  (h)  -  u>') 
( o  j  —  — 


(2tt) 


(43) 


E  (h)  =  cos  (h  a)  +  cos  (h  a)  +  cos  (h  a)  , 

A  y  Z 


(44) 


and  -  g  <  hx,  hy,  hy  <  J  .In  Fdgure  1(a)  we  plot  a  three 

hundred  histogram  approx  1  mat. Ion  to  (43)  >  whci'o  w°  have  used 


-16- 


one-mill  Ion  Monte  Carlo  points®  to  Integrate  (43).  We  notice  that 
whereas  (43)  is  completely  a  symmetric  function  around  w'  *  0,  our 
Monte  Carlo  approximation  has  a  small,  but  finite,  odd  part.  Normally 
this  imperfection  should  be  removed  by  symmetrizing,  but  here  this 
is  not  necessary  since  our  operator  form  of  the  Hilbert  kernel,  eq. 
(30),  does  not  depend  on  the  parity  of  the  expansion  functions  except 
in  that  even  functions  are  transformed  Into  odd  ones,  etc.  Moreover 
the  small  odd  part  of  g(w)  forms  a  measure  of  the  error  in  the  Monte 
Carlo  integration.  By  expanding  the  delta  function,  we  may  write 
(43)  as 


dt  cos  w* 


(45) 


a  form  first  utilized  by  Callaway.10 


The  Hilbert  transform  of  g(w 


F(w)  =  £ 


T  g(v/)  . 
\n-vi1  3 


Thus  we  have  the  analytic  representations: 


Figure  1  gives  the  result  of  integrating  the  above  forms  at  the  250 
positive  zeros  of  the  500th  order  Hermite  function  by  Hermite  integra¬ 
tion  (see  Appendix  A). 

Although  the  Hermite  expansions  are  exact  in  principle,  the 
convergence  of  the  expansion  In  Hermite  functions  may  be  substantially 
improved  by  removing  the  low  order  moments  from  g(w)  and  transforming  these 
analytically.  Then  a  Laurent  expansion  of  F(w)  will  begin  with  order 
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m+1  as 


so  that 


F(w) 


=  1  -  V  w"m  I  dw'  w' 
w  m-0  -:*> 


mg(w), 


In  practice  good  results  are  achieved  if  only  the  zeroth  moment  Is  removed, 
r'rom  g(w)  the  triangle-like  function  g  (w)  is  subtracted  where 

-  §  Q  -  SGN(w)w /£). 

This  function  is  also  plotted  in  Fig.  1  (/g$(w)dw  =  2).  The  analytic 
transform  of  this  is  given  by 


This  is  also  plotted  in  Fig.  1.  Note  that  to  first  order,  the  log  poles 
are  subtracted  out.  The  residium  of  has  been  expanded  using  the  first 
250  Hermlte  functions  and  a  scale  factor  of  a  =  where  N  »  250 

so  that,  in  effect,  the  highest  functions  just  span  the  range  from  0  to 
6.  This  scaling  avoids  any  "edge  effects"  due  to  a  confluence  of  the 
upper  level  of  the  range  of  integration  ani  the  Hermlte  function 
expansion  set.  Table  I  lists  the  first  20  expansion  coefficients. 
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Si  nee  Fig.  1  plots  both  g(u)  and  F(u)  on  the  same  scale,  we  see 
that  the  four  critical  points  of  g(w)  are  mapped  into  four  points  at 
the  same  u,  except  that  the  types1*  of  critical  points  have  inter¬ 
changed  under  Hilbert  transform;  i.e.,  M  M  M  .  We  conclude 

0  13 

from  this  example  that  critical  point  structure  cannot  shift  under 
a  proper  Hilbert  transformation. 

B.  Transition  Metal 


For  our  second  example  we  consider  the  electronic  structure  of  a 
face  centered  cubic  pure  transition  metal--here  palladium  as  calculated 
by  Mueller,  Freeman,  Diumock  and  Furdyna.  We  recall  that  the  non¬ 
interacting  electronic  Greens'  function  in  the  momentum  and  frequency 

13 

representation  for  n  bands  is  given  by 

G°(k,u)  =  E  (<o  -  E  (k)  +  i<5  SGN  (E(k)  -  y))"1  (48) 

n  n  “  n  - 

where  6  is  an  infinitesimal  and  y  is  the  Fermi  energy  Ep.  The  imaginary 
part  of  (48)  is  given  by 
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Ira  0°&.«)  ■■  -nESGN(Bn(k)  -  ^  6  (Ejk)  -  „) 


09) 


and  the  density  of  states  by: 


e(B)  -  's  SGN  («-->*)  ~r  I  im  o0(k,w)  dk 


(50) 


Tn  Figure  3(a)  we  have  plotted  the  weighted  integral  over 
i  he  Briliouin  zone  (BZ)  of  the  imaginary  part  of  G°(klU)).  The 
dispersion  relation  for  the  Greens-  function  is  given  by: 


Re  (^k>(U)  =  P  f  d(()(  Im^-,CU,)SGN(^1-^) 

17  •  -  -  P  -  (j) 


(51) 


vmcro  the  presence  of  signu,.,  function  m  (51)  removes 
the  stgnum  function  in  («),  co  that  the  integral  to  be  perfomed 
in  (51)  is  just  the  principal  part  of  a  delta  function.  Here 

W°  consider  only  the  0)  dependence  of  (51)  by  integrating  (51) 
over  tne  first  BZ  to  produce  the  Hilbert  transform  pair  G  = 


Gp  t-  iGj  where 


gr(u,)  =  f  Re  $k,w)  dk 

(Sn)0  J 


(5«) 


Gj(oj)  =  TT  g  («)  . 


(53) 


U<  have  "‘hen  Hilbert  transform  of  ('»?)  unit,,*  Ui«.  ir.*H.liocl 
i,]'  "  111  ,l  *'  i”‘l  M  l>,V  expanding  ug(u>  )  to  flr:;»  i  !  c  ,< | 

l(<  havo  ut5^  a  scale  far-..,,  a  r< .u, id 

.•ro,„W.  a  ("H  t  i. ).'■/»  wh,ro  hu„,  tor  .  -• 


notice  that  the  fact  that  g(<u)  has  no  parity  symmetry  does  not 
harm  our  completely  general  transformation,  eq.  (30).  V/e  give 
our  results  in  Figure  3(b),  and  again  list  the  first  20  expansion 
coefficients  in  Table  3.  Note  that  we  have  used  atomic  units 
of  states  per  atom-Rydberg  for  both  the  real  and  Imaginary  parts. 
The  total  time  for  our  routines  for  the  300  term  expansion  was 
1.0  minutes  on  the  IBM  360/50/75  at  the  Applied  Mathematics  Division 
at  Argonne  National  Laboratory. 

3°tisalon 

In  this  paper  we  have  considered  an  exact  numerical  procedure 
for  finding  the  Kramers-Kronlg  or  Hilbert  transform  of  a  given 
function  by  means  of  an  expansion  in  terms  of  Hermite  functions. 
Although  our  procedure  is  most  easily  applied  to  cases  where 
the  function  to  be  transformed  is  defined  over  the  whole,  real 
line,  we  have  seen  that  simple  extensions  of  1  ho  notion  of  least- 
square  fitting  allow  us  to  treat  the  proof  lea  I ly  important  case 
of  a  function  defined  only  over  a  region.  Wo  have  seen  that  the 
compactness  of  our  expansion  allows  us  to  treat;  the  special  case 
of  optical  data,  which  has  several  variables  related  by  non-linear 
equations,  by  means  of  a  single  set  of  expansion  coefficients  of 
the  Hermite  functions.  V/e  have  stressed  that  although  our  pro¬ 
cedures  are  well-defined  in  the  case  of  evenly  spaced  data,  the 
best  use  of  our  techniques  will  be  achieved  if  the  daua  is  spaced 
by  the  zeros  of  a  high  order  Hermite  polynomial,  vie  believe  that 
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1  ho  advantages  of  the  latter  procedure  are  so  great  that  all 
optical  experiments  should  be  performed  in  this  way. 

Because  of  the  simplicity  and  speed  of  expansion  of  numerical 
functions  in  terms  of  Hermite  polynomials,  we  believe  that  our 
procedures  should  prove  useful  in  a  wide  variety  of  problems. 

Moreover  we  stress  that  because  the  Hilbert 

transform  has  been  performed  analytically,  rather  than  numerically, 
our  results  are  independent  of  the  "noi3e"  or  discontinuities 
of  the  data.  Wc  believe  that  this  is  not  true  of  numerical  pro¬ 
cedures  which  operate  on  data  directly.  Finally  although  all 
of  our  discussion  of  this  paper  has  been  couched  in  terms  of  the 
Hilbert  transform,  we  point  out  that  our  procedures  are  equally 
valid  for  the  Fourier  transform  operator,  eq.  (15).  Thus  in 
contrast  to  Russel  who  concluded  in  1933  that  "this  U3e  of 
Hermite  functions  (i.o.,  Fourier  transformation)  ...  is,  in 
general,  not  practicable,",  wo  believe  that  Hermite  functions  are 
eminently  suitable  for  such  numerical  work.  Only  future  results 
can  determine  whether  such  an  expansion  technique  in  Hermite  function 
would  be  actually  faster  than  "The  Fast  Fourier  Tran,  form  Method1^." 
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APPENDIX^  E  RM ITE^  IJNT EGJl  ATI  ON 

Here  we  consider  Hermite  Integration  of  a  function  defined 
over  the  whole  real  line.  A  thorough  discussion  of  this  problem 
has  been  given  by  Ilochs Lrasser1^  and  Davis  and  Polansky^0  based 

pi  pp 

on  the  classic  work  of  Russel'  and  Salzer'  et.al.  Our  results 
differ  from  that  of  the  previous  workers  only  in  that  we  have 
gone  to  higher  order  23  and  used  a  different  weight  function.  Hermite 
integration  of  a  function  f(x)  to  order  n  is  given  by 

J  e"x  f(x)  d x  =  £  Wn(xJ)  f(x")  (Al) 

—  00  1=  1 

where  x^  is  the  1th  zero  of  the  nth  order  Hermite  polynomial  H  , 
and  the  weight  factors  W  are  given  by 


/  n. 
wn<xi> 


n; 

1,2  f'n-1 


(A2) 


-23- 


A  numerically  more  convenient  representation  for  Hermlte  inte¬ 
gration  is  achieved  if  we  consider  instead 


f.  ® 

J  e(x) 


W'n(xi>  g(xi) 


(A3) 


w lie  re  W  '  is  given  by 


W'n(x")  =  Wn  (xj1)  exp  (xj  )  =  l/(  n^_1(  xj) ) 
where  ^  is  the  n'th  order  orthonormal  Herrnite  function. 


(A*) 


In  Table  IV  we  list  the  positive  zeros  and  weight  functions  W  ' 
for  the  Hermlte  polynomials  of  order  20,  26,  50,  76,  150  and  300. 
For  the  overlapping  case  of  n  -  20,'  over  results  agree  with  those 
of  Salzer,  et;  al.  to  at  .least  14  significant  figures. 

V/e  have  found  the  acres  of  the  He  unite  polynomials  (or 
Hermlte  functions)  by  noticing  that  the  sequence  of  polynomials 
H0,  ...  Hn  form  a  Sturm  sequence.  Thus  by  counting  the  number 

of  sign  changes  In  the  sequence,  we  know  precisely,  by  Budan's 
Theorem,  the  number  of  zeros  between  a  given  point  x  and  zero, 
(infinitesimally  above  zero  the  sign  of  all  the  Hermlte  polynomial 
is  plus).  An  estimate  of  the  largest  zero  of  the  Hermlte  poly¬ 
nomial  of  order  n  is  found  by  remembering,  that  the  Hermlte 

functions  are  eigenfunctions  (w'.Mi  eigenvalue  w  =  (  m-1/2)  ftu,) 

d2  1  2 

of  the  simple  harmonic  osclll.vo  •  opn-.v  or,  H  =  !•  ^-k  x  ,  where 

w  .  we  recall  that  M  e  pro-a:  ill* density  of  :  he  nth 

harmonic  oscillator  eigenfunction  Is.  t  ou:  s<o  hy  classical 
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limit  inr  ,  or  in  the  appropriate  units  (m  =  k  =  1),  by  x  x  = 

V2n  +  1  .  Thus  all  of  the  zeros  of  the  n'th  Hermite  polynomial 

lie  between  ±  V^n  +  1. 

Our  computational  procedure  for  evaluating  the  zeros  of 
the  n'th  polynomial  used  the  Interval  given  by  the  absolute  bounds 
.fcVsn  +  1  and  continuously  sub-divided  this  interval  by  a  factor 
of  2.  Each  subinterval  was  selected  for  further  operation  depending 
on  whether  the  Sturm  sequence  predicted  an  appropriate  zero  in 
the  left  or  right-hand  portions.  Such  a  procedure  quickly,  and 
accurately  converged  to  our  results  in  Table  IV.  We  have  tested  our 
zeros  through  the  surn  rule: 

Z  K’)2  ■=  («) 

1=1  * 

to  an  accuracy  of  at  least  14  significant  figures. 


1 


\J 
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Copies  of  these  programs,  as  well  as  the  punched  cards  for 
the  zeros  are  available  on. request.  we  believe  that  because  of 
the  simple  error  spectrum  in  generating  Hermite  functions  (only 
rounding  errors),  that  our  numerical  procedures  are  accurate  to 
arbitrarily  high  order.  Cost  is  the  only  limiting  factor. 


TABLE.  I  The  first:  few  llermite  function  expansion  coefficients  of 
the  triple  cosine  density  of  states.  As  explained  in  the  text,  the 
scale  factor  a  **  6//2NI-1  v; as  used  where  here  N  l.r)0.  The  coefficient 
for  Cn  (Dn)  arc  non-zero  for  even  (odd)  order  because  they  were  de¬ 
rived  from  a  Mont e .-Carlo  function  which  was  accurate  to  ahoyt  1%. 

Note,  that  the.  convergence  of  1)  is  much  slower  than  G.  Removing  the 
zeroth  moment  from  g(w),  as  explained  in  the  text,  greatly  iiijpYover. 
this  rate  so  that  the  expansion  of  the  residual  of  1)  is  essentially 
zero  after  N  -  75. 


n 

n 

0 

1.0758 

0.0010 

1 

-0.0019 

0.2436 

2 

0.7517 

-0.0008 

3 

-0.0009 

'  0.3469 

4 

0.6029 

0.0002 

5 

-0.0042 

0.4376 

6 

0.4306 

'  -Q.0053 

7 

0.00(17 

0.4434 

•  8 

0.2842 

-0.0045 

9 

0.0037 

0.3696 

10 

0.2355 

-0.0015 

20  • 

0.1263 

(nil)  0.2815 

40 

0.0048 

0.1926 

60  . 

0.0054 

.  0.1232 

80 

-0.0002 

0.0977 

100 

0.0001 

0  .’0796 

120 

0.0008 

0.0730 

140  . 

-0.0005 

0.0693 

TABLE  IV.  The  Positive  Zeros  and  Weight  Functions  for 
Hermite  Integration  for  Various  High  Orders.  Note  that 
the  brackets  at  the  end  of  each  number  enclose  a  multi¬ 
plicative  power  of  10. 


N  =  20 

Wn<X?> 

1 

xn 

Ai 

1 

0.245340708300901 

[00' 

0.490921500666746 

(00) 

2 

0.737473723545394 

'00' 

0.493843385272053 

00) 

3 

0.123407621539532 

01 

0.499920871336291 

00) 

4 

O.I738537712II659 

01 

0.509679027117459 

00) 

5 

0.225497400208928 

01 

0.524080350943559 

00 

6 

0.27o880605842813 

01 

0.544851742364521 

00) 

7 

0.334785456738322 

01 

0.575262442852505 

00) 

8 

0.394476404011563 

01 

0.6222786961914.15 

00} 

9 

0.460368244955074 

01 

0.704332961176944 

00) 

10 

0.538748089001123 

!oij 

0.898591961453191 

(00) 

N  =  26 

wn<X?)  ■ 

i 

yn 

Ai 

1 

0.215777856243463 

00] 

O.43I682O7293I449 

[00) 

2 

0. 6460952139934481 

00 

0.433211201157933 

00) 

3 

0.1032733011077&8 

01 

O.43634.I26958I532 

00) 

4 

0.152136151665192 

01 

0.44 1226060840462 

00) 

5 

0.106585478564114 

01 

0.448,1246524005751 

00) 

6 

0.241841576477378 

01 

0.4574455005207431 

O.469827981244890I 

00) 

7 

0.288176221954309 

01 

00) 

8 

O.335942718235083 

01 

0.48629914827645H 

00 

9 

O.2S5628641990915 

01 

0.5085938683499731 

00 

10 

0.43796026o298331 

'01 

0.5399354459710701 

00) 

11 

0.494132495724138 

0.556452498195010 

01 

0. 5870633153897901 

00) 

12 

01 

0.6677573314543021 

00  ) 

13 

O.63C955036562569 

01 

0.856242029971326( 

;oo) 

N  =  50 

W,K> 

i  . 

xn 

1 

0.156302546  8  ■■■>•»  62 

CO] 

0.312630298030360 

00) 

2 

0. 46905905677' 236 

' 0 

0 . 3129334  4  >'052390 

'00) 

3 

0.  78227  .1  ( '< 1  /'  i\ ;•*'-(  >  1 

H 

0.31  •5435  6996:9642 

'00) 

4 

0.109625.11 v 

<8. 

(). ]  446  ,5(  -'2473 

’  00) 

5 

0.14113.177 >4  '  •< 

!«8.1 

0.  31572()440l 30f)')4 

’00) 

1 


6 

o. 172780654751590 

0.317316124349605 

(00) 

7 

0.204607196868641 

01 

0. 319277915993182 

00 ) 

8 

0.236649390429866 

(01 

0.321624537962057 

00) 

9 

0.268948470226775 

01 

0.324422448708867 

00) 

10 

0.301549776957452 

01 

0.327687661354607 

00) 

11 

0.334503831393789 

01 

0.331488254247523 

00) 

12 

0.367867706251527 

01 

0.335897876837687 

00) 

13 

0.401706817285813 

01 

0. 341010727258748 

00  j 

14 

0.436097316045458 

01 

0.346948769599379 

00) 

13 

0.471129366616904 

01 

0.353872469979462 

00) 

16 

O.50691175849.I724 

01 

0.361997265091619 

00) 

17 

0.543578608722495 

01 

0.3716197712499171 

00) 

.18 

0.5S129946754204K 

01 

0. 3831613921964 80' 

oo  5 

19 

0.620295251927467 

01 

0.3972449437758481 

00) 

20 

0.660864797385536 

01 

0. 4148388210590231 

00) 

21 

0.70343235097706k 

01 

0.4375532823000501 

00) 

22 

0.7486409429864201 

01 

0. 4683262119425481 

*00) 

23 

0.797562236820564( 

01 

0.5133047978515431 

00) 

24 

0.8522771030917801 

01 

0.5886052973773021 

00) 

25 

0. 9l824o6958l2932( 

I  01) 

0.7613486911180961 

'oo) 

N  =  76 

1 

x? 

1 

0.12699223100184.11 

’00] 

0.253993385676772 

(00) 

2 

0.  3'c  103026441027 6 

00 

0.25410057746354.1 

00) 

3 

0.635229306367669 

00 

0.254315550492752 

OO) 

4 

0. 869697681755121 

00 

0.254639493166670 

00 ) 

5 

0. 114 454521004 402 

01 

0.255074212122994 

00) 

6 

0.139988384 198577 

01 

0.255622161952324 

00 ) 

7 

0.165582831812922 

01 

0.256286486333430 

OO) 

8 

0.191249689274787 

01 

0.257071071989737 

00 ) 

9 

0.217001210540204 

01 

0.257980617366250 

00) 

10 

0.242850163495476 

01 

0.259020718536551 

00) 

11 

0.268809924 353835 

01 

0.260197975618256 

OO) 

12 

O.294894586891950 

01 

0.261520123961741 

OO) 

13 

0.321119078768898 

01 

0.2629961956626761 

00) 

14 

0.347499305806231 

0.1 

0. 2646367.18648907 

00) 

15 

0.374052310941034 

01 

0.2664539628697781 

00) 

16 

0.40079646366324.1 

01 

0.268:4622462074421 

00) 

17 

0.427751683211511 

(9. 

0.2706783169815881 

00) 

18 

0. 454939704768  ;C5 

01 

0.273121835  '4 76o4 1 

00) 

19 

0.4823844C  ■.  C7':  C  : 

C1 

c  .275-15  '•  '  275:  l6i 

00) 

20 

0.510112K  (V4--14C 

[Cl] 

0.27  -  Y8'-270852212( 

00) 

J 


i 

N  =  76 

Xn 

xi 

21 

0.538152412179850 

(01; 

0.282071512255670 

(00) 

22 

0.566538142750430 

01 

0.285705215988960 

00) 

23 

0.595306724459681 

01 

0.289737351776540 

00) 

24 

0.624500840679173 

01 

0.294226784433783 

00) 

25 

0.654169740880243 

01 

0.299246627607381 

00) 

26 

0.684370887971893 

01 

0.304888971325212 

00) 

27 

0.715172175595064 

01 

0.311271705500812 

0°) 

28 

0.746654979788378 

01 

0.318548640662417 

00) 

29 

0.778918468838888 

01 

0.326925003272138 

00) 

30 

0.812085877545208 

01 

0.336682064992593 

00) 

31 

0.846313978205818 

01 

0.348218086208989 

00) 

32 

0.881808021462619 

01 

0.362120212339010 

00) 

33 

0.9188466371460591 

01 

0.379299672898844 

00) 

34 

0.9578263684375551 

01 

0.40126950349675K 

00) 

35 

0.9993491992840451 

01 

0.430787054460799( 

*00) 

36 

0.1044419184013671 

01 

0.4/36l6l68889675( 

'00) 

37 

0.10949861564373H 

02] 

0.544854576829319 

00) 

38 

0.1156l5923675424( 

02] 

0.707324524548646( 

;oo) 

l 

N  -  100 

Xn 

1 

0.110795672422440 

[oo; 

0.221596255924184 

(0°) 

2 

0. 332414692342232 

0.554114823591617 

°° 

0.221650420411516 

(oo) 

3 

00 

0.221758921729024 

oo ) 

4 

0.775950741540146 

°° 

0.221922106165498 

00} 

5 

0.997977436098105 

°° 

0.222140497191774 

00) 

6 

0.122025039121895 

01 

0.222414800390545 

oo) 

7 

0.144282597021593 

01 

0.222745910173358 

00) 

8 

0.166576150-74151 

01 

0.223134918413969 

00) 

9 

O.I889II5537427O1. 

01 

0.223583125167178 

00) 

10 

O.211294799637119 

01 

0.224092051687791 

00) 

11 

0.233732046390688 

01 

0.224663456017246 

00) 

12 

O.256229640237261 

01 

0.225299351467709 

00) 

13 

0.276794142398199 

01 

0.226002028407800 

00) 

14 

O.301432356033H6 

01, 

0.226774079843686 

00) 

15 

O.32415136796310I 

01 

0. 2276l8431398403( 

00} 

16 

O.346956563641839 

01 

0. 226'538376426247( 

00 

17 

O.369661665931840 

01 

0.22953761716487K 

00) 

18 

0.392866666342767 

01 

0. 2306203 13034513( 

00) 

1$ 

0.  4159666-55  1.3103 

01 

0.23 1 79 113 7 4  53704 

00) 

20 

0. 439230r?f 

01 

0 . 23 " 055  34  4  6697 1 1 

00) 

21 

0.4626c  ■(.  ■  1' 

’01 

0.2  -44 1  501.21.7591 

00) 

22 

0.4861175C  91 79!-  8! 

[01] 

0.2  -‘5;  •  •••  322704 yoo( 

00) 

i 


23  0.509784510508914(01)  0.237471299920463(00) 

24  0.533615836013836(01)  0.23917632129956o(oO) 

25  0.557624164932992(01)  0.241013092922388(00) 

26  0.581823213520352(01  0.242992685579109(00) 

27  0.606227883261430(01)  0.245127777913569(00) 

28  0.630854436111214(01)  0.247432956126835(00) 

29  0.655720703192154(01)  0.249925086601289(00) 

30  0.680846335285880(01)  0.252623783391418(00) 

31  0.706253106o24887(oi)  0.255552000562022(00) 

32  0.731965282230454(01)  0.258736790904382(00) 

33  0.758010080785749(oi)  0.262210289443759(00) 

34  0.7844i8238446082(01)  0.266011005285956(00) 

35  0.8i1224731116279(01)  0.270185543533994(00) 

36  0.83846969404i626(01)  0.27479093832684o(oo) 

37  0.866199616813452(01)  0.279897872545095(00) 

38  0.894468921732547(01)  0.285595214531379(00) 

39  0.923342o890219i6(oi)  0.291996563960282(00) 

40  0.952896582339012(01)  0.299249958039083(00) 

41  0.983226980777797(01)  0.307552728501676(00) 

42  0.101445099412928(02)  0.317175110971854(00) 

43  0.104671854213428(02)  0.328499484645464(00) 

44  0.108022607536847(02)  0.342089262314264(00) 

45  0.111524043855851(02)  0.3588l84090711l6(oo) 

46  0.115214154007870(02)  0.380137438800590(00) 

47  0.1191506l9431l42(02)  0.408688658441919(00) 

48  0.123429642228597(02)  0.449993171054426(00) 

49  0.128237997494878{o2)  0.518506807270215(00) 

50  0.134064873381449(02)  0.674353552420927(00) 


1 

0.9053935621739301 

OH 

2 

0.2716279329434621 

00 ) 

3 

0.4527461 1652.14  30< 
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